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Abstract Effects of disorder on superconducting properties of noncentrosymmet- 
ric compounds are discussed. Elastic impurity scattering, even for scalar impurities, 
leads to a strongly anisotropic mixing of the electron states in the bands split by the 
spin-orbit coupling. We focus on the calculation of the critical temperature Tc, the 
upper critical field H^2, and the spin susceptibility Xij- It is shown that the impurity 
effects on the critical temperature are similar to those in multi-band centrosymmet- 
ric superconductors, in particular, Anderson's theorem holds for isotropic singlet 
pairing. In contrast, scalar impurities affect the spin susceptibility in the same way 
as spin-orbit impurities do in centrosymmetric superconductors. Another peculiar 
feature is that in the absence of inversion symmetry scalar disorder can mix singlet 
and triplet pairing channels. This leads to significant deviations of the upper crit- 
ical field from the predictions of the Werthamer-Helfand-Hohenberg theory in the 
conventional centrosymmetric case. 



1 Introduction 



Discovery of superconductivity in a heavy-fermion compound CePtsSi (Ref. f\\) 
has stimulated considerable interest, both experimental and theoretical, in the prop- 
erties of superconductors whose crystal lattice lacks a center of inversion. The list of 
noncentrosymmetric superconductors has been steadily growing and now includes 
dozens of materials, such as Ulr (Ref. g]), CeRhSis (Ref. 0), CeIrSi3 (Ref. H), 
Y2C3 (Ref. Ill), Li2(Pdi_.,^,Pt,03B (Ref. 11), and many others. 

A peculiar property of noncentrosymmetric crystals is that the spin-orbit (SO) 
coupling of electrons with the crystal lattice qualitatively changes the nature of 
the Bloch states, lifting the spin degeneracy of the electron bands almost every- 
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where in the Brillouin zone. The resuhing nondegenerate bands are character- 
ized by a complex spin texture and a nontrivial wavefunction topology in the 
momentum space. Q This has profound consequences for superconductivity, in- 
cluding unusual nonuniform superconducting phases, both with and without mag- 
netic field,|l8l m [lOl HU [El [13 magnetoelectric effect, Ql US [H HT] [TSl and 
a strongly anisotropic spin susceptibility with a large residual component at zero 
temperature. ifTSl [191 l20l 1211 These and other properties are discussed in other chap- 
ters of this volume. 

In this chapter we present a theoretical review of the effects of nonmagnetic 
impurities in superconductors without inversion symmetry. In Sect. [2] the disorder- 
averaged Green's functions in the normal and superconducting states are calculated. 
In Sect. [3] the equations for the superconducting gap functions renormalized by 
impurities are used to find the critical temperature T^. In Sect. [H the upper critical 
field Hi-2 is calculated at arbitrary temperatures. In Sect. [Sj we calculate the spin 
susceptibility, focusing, in particular, on the effects of impurities on the residual 
susceptibility at T = 0. Sect.|6]contains a discussion of our results. Throughtout this 
chapter we use the units in which kg = h= \ . 



2 Impurity scattering in normal and superconducting state 

Let us consider one spin-degenerate band with the dispersion given by eo(^), and 
turn on the SO coupling. The Hamiltonian of noninteracting electrons in the pres- 
ence of scalar impurities can be written in the form H ^ Hq+ Hi„,p, where 



a , j3 =t, I is the spin projection on the z-axis, Y,k stands for the summation over the 
first Brillouin zone, a are the Pauli matrices, and the chemical potential is included 
in eo(^). The "bare" band dispersion satisfies eo(— ^) = eo(^), £o{8^^k) = £o{k), 
where g is any operation from the point group G of the crystal. The electron-lattice 
SO coupling is described by the pseudovector y(ife), which has the following sym- 
metry properties: y(A:) = — y(— it), gy(g^'^) = Y{k)- Its momentum dependence 
crucially depends on G, see Ref. [T]- For example, in the case of a tetragonal point 
group G = C4,,, which is realized, e.g., in CePtsSi, CeRhSis, and CelrSis, the sim- 
plest expression for the SO coupling is y(ik) = Yo{kyX — k^y), which is also known as 
the Rashba model. [22] In contrast, in a cubic crystal with G = O, which describes 
Li2(Pdi_;„Pt^)3B, we have y(*) = yok. 

The impurity scattering is described by the following Hamiltonian: 



(1) 
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The impurity potential U{r) is a random function with zero mean and the correlator 
{U {ri)U {r2))imp = nimpUQ5{r\ — r2), where nimp is the impurity concentration, and 
f/o is the strength of an individual point-like impurity. 

The Hamiltonian ([T]) can be diagonalized by a unitary transformation ata — 
Ex "a/l {^)ckX^ where /I = ± is the band index (helicity), and 

with the following result: 

it A=± 

The energy of the fermionic quasiparticles in the Ath band is given by t,x (k) ~ 
Eo{k) +X\Y{k)\. This expression is even in k despite the antisymmetry of the SO 
coupling, which is a manifestation of the Kramers degeneracy: the states \kX) and 
I — ^, A) are related by time reversal and therefore have the same energy. In real non- 
centrosymmetric materials, the SO splitting between the helicity bands is strongly 
anisotropic. Its magnitude can be characterized by Esq = 2 max^t I y(^) I ■ For in- 
stance, in CePtaSi Eso ranges from 50 to 200 meV (Ref. 1231 ). while in Li2Pd3B it 
is 30 meV, reaching 200 meV in Li2Pt3B (Ref. Il24l ). 

In the band representation the impurity Hamiltonian (|2]l becomes 

where f is the system volume, U{q) is the Fourier transform of the impurity po- 
tential, and w{k,k') = {f{k)u{1^). We see that the impurity scattering amplitude 
in the band representation is momentum-dependent, even for isotropic scalar im- 
purities, and also acquires both intraband and interband components, the latter 
causing mixing of the helicity bands. In the case of a slowly-varying random po- 
tential, keeping only the forward-scattering contribution U [q) ^ dq o, one obtains: 
wxxi{k,k') ~ dxx'^k i-^- bands are not mixed. 

The electron Green's function in the helicity band representation is introduced 
in the standard fashion: Gxx'i^^ x;k' = — {TT:Ckx{'^)c\i xX'^')) ■ in ti^s absence of 
impurities, we have G^^ xX'i^i C^n) = 8xx'/[i^n — ^xi^)]^ where (0„ = (2« + \)kT is 
the fermionic Matsubara frequency. 

We will now show that the impurity-averaged Green's function remains band- 
diagonal. The disorder averaging with the Hamiltonian (|5]l can be performed using 
the standard methods, ||25ll resulting in the Dyson equation of the form G ' = Gq ' — 
£, where G is the average Green's function and I. is the impurity self-energy, see 
FigH in the Born approximation, taking the thermodynamic limit ^ ^ oo, we have 

t{k,(o„) = «,-,„Xy j^w{k,k!)G{k! ,(On)w{k ,k). (6) 



K. V. Samokhin 



/ 



Fig. 1 The impurity self-energy in the band representation. The dashed line corresponds to ni,„pUQ, 
the vertices include the anisotropy factors w{k, k'), and the solid line is the average Green's function 
of electrons in the normal state. It is shown in the text that the self-energy is nonzero only if Xi = A4 
and A2 = A3. 



Seeking solution of the Dyson equation in a band-diagonal form: Gui = Gx dxx'i 
the integrand on the right-hand side of Eq. (|6]) can be written as follows: 

where f,- are the Pauli matrices, and y = Y/\Y\- The second line in this expression 
vanishes after the momentum integration, therefore E{k,co„) ~ Z(a)„)fo. The real 
part of the self-energy renormalizes the chemical potential, while for the imagi- 
nary part we obtain: lmE{(0„) = — Fsigno,,. Here F = miimpUQNf is the elas- 
tic scattering rate, with A^^ defined as follows: Np — {N+ +N-)/2, where A^;^ ~ 
' Y,k 5 [^i {k)] is the Fermi-level density of states in the Ath band. Thus we arrive 
at the following expression for the average Green's function of the band electrons: 

Gxx'ik,(o„)= . . (7) 

This derivation is valid under the assumption that the elastic scattering rate is small 
compared with the Fermi energy which justifies neglecting the diagrams with 
crossed impurity lines in the self-energy in Fig.[T] 



2.1 Impurity averaging in superconducting state 

In the limit of strong SO coupling, i.e. when the band splitting exceeds all super- 
conducting energy scales, the Cooper pairing between the electrons with opposite 
momenta occurs only if they are from the same nondegenerate band. The pairing 
interaction in the strong SO coupUng case is most naturally introduced using the 
basis of the exact band states. |fT9l l23l l26l which already incorporate the effects of 
the crystal lattice potential and the SO coupling. The total Hamiltonian including 
the pairing interaction is given hy H = Hq+ Hj„,p + //,„,, where the first two terms 
are given by Eqs. ([T]i and Q respectively, and the last term has the following form: 
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'^-k' .k'^k'+q.X'- 



(8) 



kk'qXX' 



Physically, the pairing interaction is mediated by some bosonic excitations, e.g. 
phonons, and is effective only at frequencies smaller than a cutoff frequency (Oc, 
which has to be included in the appropriate Matsubara sums. Alternatively, the cut- 
off can be imposed on the momenta in Eq. ([8]), as in the original Bardeen-Cooper- 
Schrieffer (BCS) model. The diagonal elements of the pairing potential Vxx' de- 
scribe the intraband Cooper pairing, while the off-diagonal ones correspond to the 
pair scattering from one band to the other 

The pairing potential can be represented in the following form: Vxxi{k,k') — 
tx{k)tli{k)Vxx'{kik')^ see Ref. |j27l. Here tx{k) = -tx{-k) are non-trivial phase 
factors originating in the expression for the time reversal operation for electrons 
in the helicity bands: K\kX) = tx{k)\ —k,X), 1191 l26l while the components of 
Vxx' are even in both k and k' and invariant under the point group operations: 
Vx}ii{g~^k,g^^k') = Vxx'ik,k'). The latter can be expressed in terms of the basis 
functions of the irreducible representations of the point group. [281 In general, the 
basis functions are different for each matrix element. Neglecting this complication, 
and also considering only the one-dimensional representation corresponding to the 
pairing channel with the maximum critical temperature, one can write 



where the coupUng constants Vxx' form a symmetric positive-definite 2x2 ma- 
trix, and ^xi^^) ^re even basis functions. While 0+(fc) and have the same 
symmetry, their momentum dependence does not have to be the same. The ba- 
sis functions are assumed to be real and normalized: = 1' where the 
angular brackets denote the Fermi-surface averaging in the Ath band: {{■■■))x = 



Treating the pairing interaction ^ in the mean-field approximation, one intro- 
duces the superconducting order parameters in the helicity bands, which have the 
following form: Ax {k, q) = tx ik)Ax {k,q). The superconducting order parameter is 
given by a set of complex gap functions, one for each band, which are coupled due 
to the interband scattering of the Cooper pairs and other mechanisms, e.g. impu- 
rity scattering. Thus the overall structure of the theory resembles that of multi-band 
superconductors. II29I l30l If the pairing corresponds to a one-dimensional represen- 
tation, see Eq. (|9]l, then we have Ax {k, q) = rjx iq)<j>x (*)■ 

Important particular case is a BCS-like model in which the pairing interaction is 
local in real space: 



Vxx'{k,k') = -Vxx'MkWv{k'), 



(9) 



{i/Nxr)Lk{-miik)]- 




(10) 



where y > is the coupling constant. One can show 1271 that in this model there is no 
interband pairing for any strength of the SO coupling, the order parameter has only 
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one component 77, the gap symmetry corresponds to the unity representation with 
(ik) = 1, and all coupling constants in Eq. Q take the same value: YxX' = 
Let us calculate the impurity-averaged Green's functions in the superconducting 
state. To make notations compact, the normal and anomalous Green's functions ll25l 
can be combined into a 4 x 4 matrix '^(k\,k2\'i) = —iJzCk^(x)C\^^(^)), where 

Ck = (cjfcAiclj^;)^)^ are four-component Nambu operators. Averaging with respect 
to the impurity positions restores translational invariance: {'^^{ki,k2',C0„))i„,p = 

and the hats denote 2x2 matrices in the band space. The average matrix Green's 
function satisfies the Gor'kov equations: — Ei,„p)^ = 1, where 

^o-V.,co,0=f%f(f^-f?) V (12) 

and the impurity self-energy in the self-consistent Born approximation is 

f d^k' 

E„„p{k,co„) = ni,„pUiJ ■j^Wik,k')'^{k\(On)Wik\k), (13) 

which is the Nambu-matrix generalization of Eq. The 4x4 matrix W is defined 
as follows: W{k,k') = diag [v!>(^,^'), — m>-^(— —k)]. It is straightforward to show 
that [w'^{—k', —k)]xx' = f| {k)t\i{k')wix'{l^,l^')- We assume the disorder to be suffi- 
ciently weak, so that it is legitimate to use the Born approximation. Although there 
are some interesting qualitative effects in the opposite limit of strong disorder, such 
as the impurity resonance states. ||3T1 these are beyond the scope of our study. 

In the absence of impurities, the Green's functions have the following form: 
Go,Av(*::««) = hx'Go,),{k,a)n), andpQxx'ik^cOn) = 5iiiti{k)Fo^i{k,a)n), where 

In the presence of impurities, we seek solution of the Gor'kov equations in a band- 
diagonal form and require, for consistency, that the Nambu matrix components of 
the self-energy are also band-diagonal. Then, 

/Z" (*:,a)„)i:||, \^ / Ly{(o„) tx{k)i:2{(0„) 
\El\,{k,(o„) Zf^,{K(On) ) \tl{k)mcOn) -Er{-co„) 

where Ei and E2 satisfy the equations 



Effects of impurities in noncentrosymmetric superconductors 7 

1 f (fik 

1 f d^k ^ ' 

M(On) = :^n„npUiY, J (^^^(*^'®«)- 

Absorbing the real part of Ei into the chemical potential, we have Ei {co„) = iEi (©„), 
where Ei is odd in (0„. 

Solving the Gor'kov equations we obtain the following expressions for the 
disorder-averaged Green's functions: 

^ /t „ A /d)„ + ^;t(A:) 



d),?+^l(*)+li)A(*,««)p' 

(17) 



where (b„ = (0„ — Ei (o),,) and Dx {k, (O,,) = Ax {k) + E2{(0„). Substituting these into 
Eqs. (fT6b . we arrive at the self-consistency equations for the Matsubara frequency 
and the gap functions renormalized by impurities: 

W„ = ft)„ + ^^p;t(— ^== 1^" ) , (18) 

where 

p± = -^ = \±5 (20) 

are the fractional densities of states in the helicity bands. The parameter 5 = (A^+ — 
N-)/{N+ +N-) characterizes the strength of the SO coupling. 



3 Gap equations and the critical temperature 

The Gor'kov equations must be supplemented by self-consistency equations for the 
order parameter components, which have the form usual for two-band supercon- 
ductors. In particular, for a uniform order parameter, rix{q) = rjxS{q), we have 
^/L' ^il'^X' — TY,„ ik^X (^7 <J^n)<j>x (^) (recall that the basis functions are assumed to 
be real). Using Eq. ilT\ . we obtain: 

y:vxi,m'=^PxNFTy'(^^t^^t^) . (21) 
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These equations are called the gap equations and, together with Eqs. ( fTST l and (fT9] l. 
completely determine the properties of disordered noncentrosymmetric supercon- 
ductors in the uniform state. The prime in the Matsubara sum means that the sum- 
mation is limited to co„ < (Oc, where cOc is the BCS frequency cutoff. 

The superconducting critical temperature can be found from Eq. jTH after lin- 
earization with respect to the order parameter components. It follows from Eq. 
( fTSl ) that a)„ — CO,, +Fsigna)„ near Tc, and we obtain from Eq. (l% that Z2(fi)„) = 
(F /2 1 0),, I ) P A ( ) ^ A (here and below we omit, for brevity, the arguments of the 
basis functions and the subscripts A in the Fermi-surface averages). Therefore the 
linearized gap equations take the form Y,x' ^XX'^X' = 0> where 

1 r 

axx' = j^yxli-Px^xx'SiM - -^PxPx'{h){'l>x')Sn, (22) 

with Ski — ^TY,n \(')n\^'' {\(0„ \ + \r)^' . The Matsubara sums here can be easily cal- 
culated: 

^ co„+r kT \r 

where C ~ 0.577 is Euler's constant, 

+ (23) 

and f (x) is the digamma function. Note that the expression ( l23T i for the impurity 
correction to 5oi is valid if F ^ (Oc , when it is legitimate to extend the summa- 
tion in 2nTY^„[\ / {(On +F) — I /(£)„] to infinity and express the result in terms of 
the digamma functions, see Ref. |32J. Similarly, we obtain: = ^{T jT^jT . It 
is convenient to introduce the following notation for dimensionless coupling con- 
stants: 

&XX' ^ NfVxx'Px' = yxx'Nx' (24) 

(note that the matrix g is not symmetric, in general). Then, the superconducting 
critical temperature 71 is found from the equation gM) ~ 0, where 

Mxx' = -8xx' In ^ + {8,x' ^ ('^;i> (-^A')) ^ (j) ' ^^^^ 



see Ref. | 

In the absence of impurities, the second term in M^x' vanishes, and we obtain 
the critical temperature of a clean superconductor: 

Tci) = e '^ (26) 

7Z 

where 

(27) 
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is the effective coupling constant. In the presence of impurities, the cases of conven- 
tional and unconventional pairing have to be considered separately. 

Unconventional pairing. In this case {<^x) = 0, and we obtain the following equa- 
tion for Tc'. 

lnZ£^ = ^f5y (28) 



Tc \r 

The reduction of the critical temperature is described by a universal function, which 
has the same form as in isotropic centrosymmetric superconductors with magnetic 
impurities, f34| or in anisotropically paired centrosymmetric superconductors with 
nonmagnetic impurities. 13511281 In particular, at weak disorder, i.e. in the limit F <C 
Tco, we have Tc = Tco — TtF /4. The superconductivity is completely suppressed at 
U = (7r/2eC)r,o. 

Conventional pairing. Assuming a completely isotropic pairing with = 1, we 
obtain: 

InZ^. \ (29) 

Tc C2 + C3^(x) + ^yC4+C5^{x)+C6^'^{x) g 

where x^Tc/F, and 

P+{g—-g+-) + P-{g++-g-+) g+++g— det| 

Cl= ^ , C2 = , C3==— , 

C4 = ^ '^^^ 2'^" ) +^+-^- + ' C5 = (C2-Cl)det|, C6 = C3. 

We see that the critical temperature depends on nonmagnetic disorder, but in con- 
trast to the unconventional case, the effect is not described by a universal Abrikosov- 
Gor'kov function. 1361 [33I At weak disorder the suppression is linear in the scatter- 
ing rate, but with a non-universal slope: 



Tc^TcQ-- 



K 1 / , C5 
C1-7- C3 + -— 



(30) 



In the case of strong impurity scattering, F Tco, we use ^(x) = ln(l /x) + 0(1) 
at jc — > 0, to find that the critical temperature approaches the limiting value given by 

r; = r,oexpQ-^^, (3i) 

i.e. superconductivity is not completely destroyed by impurities. The explanation 
is the same as in the conventional two-gap superconductors, see e.g. Refs. I37ll38l 
|39l : Interband impurity scattering tends to reduce the difference between the gap 
magnitudes in the two bands, which costs energy and thus suppresses Tc, but only 
until both gaps become equal. One can show that both the coefficient in front of F 
in Eq. ( l30l l and the exponent in Eq. ( |3T1 ) are negative, i.e. T* < Tcq. 
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In the BCS-like model ( fTOl i the pairing is isotropic and described by a single 
coupling constant Y^X' = ^/2, and we have g = NpY . Although the expression ( [26] l 
for the critical temperature in the clean case has the usual BCS form, the analogy 
is not complete, because the order parameter resides in two nondegenerate bands, 
and Np = (A^+ +N-)/2. In the presence of impurities, the right-hand side of Eq. 
(I29] ) vanishes, therefore there is an analog of Anderson's theorem: The nonmagnetic 
disorder has no effect on the critical temperature. 

Another important particular case, possibly relevant to CePtsSi, is the model in 
which only one, say X = +, band is superconducting, while the other band remains 

normal. 1231 [33)1 This can be described by setting V+- =V =0. Using Eq. (IZSl) . 

we obtain the following equation for the critical temperature: 



where cq = 1 — p+{(j)+)^/2. At weak disorder we have Tc = Tcq — co(7rF/4), while 
at strong disorder Tc = TcQ{nTcQ / e*^ FY I '^^ . If the pairing is anisotropic but conven- 
tional, then, unlike the unconventional case with (0+) ~ 0, the superconductivity is 
never completely destroyed, even at strong disorder. 



3.1 Isotropic model 

Finding the superconducting gap at arbitrary temperatures and impurity concen- 
trations from the nonlinear gap equations ( 1211 1 is more difficult than the calcu- 
lation of Tc- We focus on the case when the pairing is completely isotropic, i.e. 

ik) = ([)- {k) = 1 and Ai {k) = rji . The order parameter components can be cho- 
sen to be real, and the gap equations take the following form: 



We further assume that the difference between p+ and p_ can be neglected and the 

pairing strength, see Eq. (|9]), does not vary between the bands: V++ = V > 0. For 

the interband coupling constants, we have = ^-+- The gap equations have two 

solutions: 17+ = i]- = 17 and 77+ = — Tj- = rj.ln the spin representation, the former 
corresponds to the singlet state, while the latter - to the "protected" triplet state. 1201 
The impurity responses of these two states turn out to be very different. 

T7+ = T7_ = 77. In this case D+(co„) — D_{(On) = D{(On), and Eqs. ( fTSl l and ( fT9] l 
take the following form: 




(32) 




(33) 



D 
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Introducing Z{cOn) = 1 + F / co]; + 1]'^ , the solution of these equations is D{cOn) — 
Z{(On)ri, Shi = Z{(0„)(0n. Therefore, the gap equation ( [33] l becomes 

where g\ = + V+-)Nf- The scattering rate has dropped out, therefore there is 
an analog of Anderson's theorem: neither the gap magnitude nor the critical temper- 
ature, are affected by impurities. Namely, we have Tc {r) = T^-q, see Eq. ( |26] | with 
g ~ gi, while the gap magnitude at T = is given by the clean BCS expression: 
r7(r=0) = r7o = (7r/eC)r,o. 

TJ+ = — Tj_ = 77. In this case D+{cOn) = —D^{(On) ~ r\, and we obtain fromEqs. 
(1181), (1191), and ([331): 

d)„ = ft)„+r^=^L=, ^^Kg2TY^'^^==, (35) 

where g2 = — V+-)Nf. In the absence of impurities, the critical temperature 
is given by the BCS expression ( |26] l with g — gi- If V+- > (attractive interband 
interaction), then g2 < g\ and the phase transition occurs into the state T7+ = Tj-. 
However, if V+_ < (repulsive interband interaction), then g2 > gi and the phase 
transition occurs into the state 77+ = — Tj_. In contrast to the previous case, both 
the critical temperature and the gap magnitude are now suppressed by disorder The 
former is determined by the equation (|29] l, which takes the same universal form as 
the Abrikosov-Gor'kov equation ( |28] l. Superconductivity is completely destroyed if 
the disorder strength exceeds the critical value = (;r/2e'')rco. 

To find the gap magnitude at T = as a function of F we follow the procedure 
described in Ref. |i40l. Replacing the Matsubara sum by a frequency integral in the 
second of equations ( 1351 ), we obtain: 




1 



2 



(36) 



where 770 = 217 is the B CS gap m agnitude in the clean case, and &> satisfies the 
equation d) = © + Ffij/ \/ &P- + rj^. Transforming the second term on the right-hand 
side of Eq. ( |36] | into an integral over a) we arrive at the following equation: 



ln-p = — -ln(2x) + 0(x-i; 



- - arctan \/ — I — — - 

2 2x 



ln{x+ \/x^-l) 

(37) 



where x = F /rj. This equation does not have solutions at F > i7 , which is consis- 
tent with the complete suppression of superconductivity above the critical disorder 
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Strength. In the weak disorder limit, F ^ f^, the solution is ~ r/lFc, while at 
r ^ I^. we have x ~ ^rc/12(r,. - F). 



4 Upper critical field at arbitrary temperature 

In this section, we calculate the upper critical field HdiT) of a disordered noncen- 
trosymmetric superconductor described by the BCS-like model ( fTOb . We assume a 
uniform external field H and neglect the paramagnetic pair breaking. The noninter- 
acting part of the Hamiltonian is given by 

h^£^{K) + y{K)a + U{r), (38) 

where K — —/V + (e/c)A(r), A is the vector potential, and e is the absolute value 
of the electron charge. The superconducting order parameter in the model (fTOb is 
represented by a single complex function r\ (r). According to Sect. |3] the zero-field 
critical temperature is not affected by scalar impurities. The critical temperature at 
H y^Q,OT inversely the upper critical field as a function of temperature, can be found 
from the condition that the linearized gap equation [V^' — TY,„^X{(Oi,)]rj{r) = has 
a nontrivial solution. Here the operator X(co„) is defined by the kernel 

X{r,r';co„) = i(tr|+G(r,r';a)„)|G^(r,r';-to„))„„^^, (39) 

where g — 102- The angular brackets denote the impurity averaging, and G(r,r'; ©„) 
is the Matsubara Green's functions of electrons in the normal state, which satisfies 
the equation (/©„ — li)G{r,r'; co„) = S{r— r'), with h given by expression ( l38b . 

The impurity average of the product of two Green's functions in Eq. (|39] l can be 
represented graphically by the ladder diagrams, see Fig.|2](as before, we assume the 
disorder to be sufficiently weak for the diagrams with crossed impurity lines to be 
negligible). In order to sum the diagrams, we introduce an impurity-renormalized 
gap function D(r, O),,), which a matrix in the spin space satisfying the following 
integral equation: 

D{r,co„)^riir)g 

+ ^riimpU^g J cPr' tig^ {G{r, r'; (0„))i,„pD{r' , (0„ ) (G^ (r, r'; - «„)) imp 

+ ^ni,„pUig J d^r' tig^G{ry ■,(o„))i,npb{r' ,(o„){G'^ {r,r' ■^(On))i,np. (40) 

This can be easily derived from the ladder diagrams in Fig.|2] by representing each 
"rung" of the ladder as a sum of spin-singlet and spin-triplet terms: 

riimpUoS^vSpa = ^riimpUog^ipgly + ^«™pt/o5fipf (41) 
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Fig. 2 Impurity ladder diagrams in the Cooper channel. Lines with arrows correspond to the av- 
erage Green's functions of electrons, g = idi, and the impurity (dashed) line is defined in the text, 
see Eq. <4U . 



where g = i&G2- 

Seeking solution of Eq. (|40] | in the form D(r, O),,) = dQ{r, (0„)g + d{r, (0„)g, we 
obtain a system of four integral equations for da, where a ~ 0, 1 , 2, 3; 

3 

^ [dab - r^abicon)] db{r, «„) = 7] (r)5„o. (42) 

b=0 

Here the operators ^^f,(a)„) are defined by the kernels 

Xb{r,r';(0„) = j^tTgl{Gir,r';co„))i,„pgb{G^iry;-co„))„„p, (43) 

with go = g, and g„ = ga for a = 1,2,3. We see that, in addition to the spin-singlet 
component do (r, a)„ ), impurity scattering also induces a nonzero spin-triplet compo- 
nent d{r, (On). The linearized gap equation contains only the former. Indeed, using 
Eqs. ( |42] | we obtain: 

W - ^tT M'-^''"^''^^''^ = 0. (44) 

NfV „ r ^ ' 

It is easy to see that the triplet component does not appear in the centrosymmetric 
case. Indeed, in the absence of the Zeeman interaction the spin structure of the 
Green's function is trivial: G„^(r, /■';©„) = 5apG{r,r' \(0n)- Then it follows from 
Eq. (US that %h{oin) = 5ab^((0„), therefore do = {l-Tfy^ri and d = 0. 

The next step is to find the spectrum of the operators ^^^(a),,). At zero field, the 
average Green's function has the following form: 

Goik, (On) = £ i7;i (fc)G;i {k, «„), (45) 

where 17^ = (1 + Ayd')/2 are the helicity band projection operators, and G^ {k, (0„) 
are the impurity-averaged Green's functions in the band representation, see Eq. (|7]). 
At H ^ 0, we have {G{r,r';(On))imp = Go(r-r';a)„)exp[(ie/c) j'^ Adr], where the 
integration is performed along a straight line connecting r and r M4l1 This approxi- 
mation is legitimate if the temperature is not very low, so that the Landau level quan- 
tization can be neglected. It follows from Eq. (03]) that ^^,i(a)„) = Yab{q, ('^n)\^^jj, 
where D = —iV + (2e/c)A and 
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Yabiq,0}„) = J j^tTglGQ{k + q,co„)gbGli-k,-o}„). (46) 

Substituting here expressions ( |45] ) and calculating the spin traces, we obtain for the 
singlet-singlet and singlet-triplet contributions: 

Y^=2^P>^{\co„\+r + iv^ {k)q sign co„/2 ) ' ^^^^ 

Yo, = i-zo = ^ E Ap 7 , . ^ , .^'^''l . -r y (48) 

2Y \|a)„|+r + iV;L(*)«signa)„/2/ 

where Vx = d£,x/dk is the quasiparticle velocity in the Ath band. We see that the 
singlet-triplet mixing occurs due to the SO coupling and vanishes when p+ = p_ = 1 
and v+ = v_ = Vf. The triplet-triplet contributions can be represented as follows: 

(1) _ 1 / l^-WfiW 



yf'' = -Vo / riy^ujy-^) \ ..g. 

'j 2j;^^^\\(0„\+r + ivxik)qsignco„/2/' ^ ' 



and 



The singlet impurity scattering channel, which is described by the first term in Eq. 
(|4TI ). causes only the scattering of intraband pairs between the bands. In contrast, 
the triplet impurity scattering can also create interband pairs, which are described 

byy,f. 

It is easy to show that if the SO band splitting exceeds both (Oc and F, then the 
interband term in Yij is smaller than the intraband one. Let us consider, for example, 
isotropic bands with ^±{k) = eo(fc) ± Y- Neglecting for simplicity the differences 
between the densities of states and the Fermi velocities in the two bands and setting 
^ = 0, we obtain from Eqs. ( |49] | and (BOl i: 



^i''^^''^") - 3(|a)„|+r)(l + .^) " ^"-(«")^- 

where r(a)„) = yj (|fi)„| Due to the BCS cutoff, the maximum value of ©„ in 
the Cooper ladder is equal to o); , therefore r„„„ ^ E^o I max(a)c, F). We assume that 
this ratio is large, which is a good assumption for real materials, therefore 



Y inter ( ) 

max - 



2 



max(tOc,r)"' ^ 



«1, (51) 
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at all Matsubara frequencies satisfying \(o„ \ < cOc- 

Thus the interband contributions to the Cooper ladder can be neglected, and we 
obtain: 

4',c^ = tLPa ( 77-7-7^-— 777^ — ■■ — 7-1^), (52) 
2Y \|fi)„|+r + ;v;L(A:)^signa)„/2/ 

where Ai Q{k) = 1 and Ai ^ik) = X%{k) for a = 1,2,3. Making the substitution 
q —T' D,we represent as a differential operator of infinite order; 

€bi(0„) = i £du £p;l (A;,„(*:)AA,i,W^^-""'^'*'''"^"'""/') . (53) 

In order to solve Eqs. ( |42] |. with the operators ^%fo(e)„) given by expressions 
$53i . we follow the procedure described in Ref. B2II . Choosing the z-axis along 
the external field; H = Hz, we introduce the operators a± = ^//(£)v± iDy)/2, and 
03 — £hD^, where £h = ^ cjeH is the magnetic length. It is easy to check that 
a+ = fl^ and [a_ , 0+] = 1 , therefore a± have the meaning of the raising and lowering 
operators, while 03 = aj commutes with both of them; [fl3,fl±] = 0. It is convenient 
to expand both the order parameter 77 and the impurity-renormalized gap functions 
da in the basis of Landau levels \N ^p), which satisfy a+|A^,/:)) = \/N + 1\N + l,p), 
a-\N,p) = \/N\N-1,p), and a3\N,p) = p\N,p), where = 0, 1, and pis a real 
number characterizing the variation of the order parameter along the field. We have 

Tl{r) = Y,riN.p{r\N,p), da{r,C0n) ^Y.''N.pi^n){r\N,p). (54) 

N.p N.p 

According to Eqs. ( l42l i. the expansion coefficients can be found from the following 
algebraic equations; 

5ab5m'Spp'-J^{N,p\?K,hi(^)\N',p') d'^,p,{(o„)^5MriN,,,. (55) 

N'.p'.b 

Substituting the solutions of these equations into 

— 77^,p - KTj^ = 0, (56) 

see Eq. (l44l i. and setting the determinant of the resulting linear equations for rj^ p to 
zero, one arrives at an equation for the upper critical field. 



4.1 Hc2{T) in a cubic crystal 



In the general case, i.e. for arbitrary crystal symmetry and electronic band structure, 
the procedure outUned above does not yield an equation for H^iiT) in a closed 
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form, since all the Landau levels are coupled. In order to make progress, we focus 
on the cubic case, G = O, with a parabolic band and the SO coupling given by 
y{k) = Yok. As for the parameter 5, which characterizes the difference between the 
band densities of states, see Eq. ( l20b . we assume that 



5c<|5|<l, 



(57) 



where 8^ = max{cOc,r)/ef <C 1. While the first inequality follows from the condi- 
tion ( fSTl i. which ensures the smallness of the interband contribution to the Cooper 
impurity ladder, the second one is always satisfied, with 1 5 1 ^ 1 corresponding to 
the rather unrealistic limit of extremely strong SO coupling. 

In order to solve the gap equations, we make a change of variables in the triplet 
component: d± ~ {di ± id2)/\/2. Then, Eqs. (|42] | take the following form: 



-r%3 i-r#33 



-r#3_ 



-r#3+ 



\ 



















d+ 







/ 


V-J 







(58) 



where 



■o± 



'02 



V2 



3f+ 



(59) 



with -S^ab = %a given by Eqs. 

According to Eq. dSSl l, one has to know the matrix elements of the operators 
?Ki,{(On) in the basis of the Landau levels \N,p). After some straightforward algebra, 
see Ref. 1431 for details, we find that ?Pqo, ^%3, ^^3, and are diagonal in the Lan- 
dau levels, while for the rest of the operators ( |59] l the only nonzero matrix elements 
are as follows: {N , p\%^\N + 1 , p) = {N + l,p\%+\N,p), {N , p\m,^\N + I , p) = 
{N+ l,p\^3+\N,p}, and {N,p\^-\N + 2,p) = {N + 2,p\S'+\N,p}. Therefore, the 
Landau levels are decoupled, and for 77 (r) = 77 {r\N, p) (77 is a constant) the solution 
of Eqs. dSSl ) has the following form: 



I do{r,co„) \ 
di{r,co„) 
d+{r,(On) 

\d-{r,a)n) J 



/ d'^^{(0„){r\N,p) \ 
dlJco„){r\N,p) 
d^j„{co„){r\N+l,p) 
\d^^,,{(0„){r\N-l,p)J 



(60) 



At arbitrary magnitude of the SO band splitting, the singlet-triplet mixing makes 
the equation foiHc2{T) in noncentrosymmetric superconductors considerably more 
cumbersome than in the Helfand-Werthamer problem, 021 even in our "minimal" 
isotropic model. It is even possible that, at some values of the parameters, the max- 
imum critical field is achieved for N > and p ^0, the latter corresponding to a 
disorder- induced modulation of the order parameter along the applied field. Leaving 
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these exotic possibilities aside, here we just assume that N ~ p = 0. Then the only 
nonzero components of the impurity-renormalized gap function, see Eq. ( |54l i. are 
4,0 and d+Q. 

It is convenient to introduce the reduced temperature, magnetic field, and disorder 
strength: 

t-— h-— — 

Tco ' Ho ' kTm ' 

where Hq = <t>o/n^Q, <Pq ~ nc/e is the magnetic flux quantum, and = Vf/2nTco 
is the superconducting coherence length (vp is the Fermi velocity). Using the ex- 
pression for the critical temperature to eliminate both the frequency cutoff and the 
coupling constant from Eq. (|56] |. we arrive at the following equation for the upper 
critical field hdit): 



In- 



1 



2E 

;i>0 



1 



2n + l (1-Cw„)(1-Cp„) + C252^2 



(61) 



where 



Pn = 



qn 



Jo 

/■I 1 - t2 r 
dpe-^"P / ds-— 
Jo 2 



2(l-.2)/4 



(62) 



dpe-'"P ( dsJ-p{l-s')e-'P'^'-^'y\ 
Jo V 4 



where 9„ = {2n + l)t+^. 

In the clean limit, i.e. at 0, or if the SO band splitting is negligibly 
small, i.e. at 5 — )- 0, the Helfand-Werthamer equation for Hc2 in a centrosymmetric 
superconductor ll42l is recovered. Therefore, the absence of inversion symmetry af- 
fects the upper critical field only if disorder is present. One can expect that the effect 
will be most pronounced in the "dirty" limit, in which Eq. i6T[ takes the form of a 
universal equation describing the magnetic pair breaking in superconductors: 1441 



(63) 



Here the parameter 



2 + 52 



(64) 



characterizes the pair-breaker strength. Note that the corresponding expression in 
the centrosymmetric case is different: Gcs — Analytical expressions for the 

upper critical field can be obtained in the weak-field limit near the critical tempera- 
ture: 

24 C 



{2 + 5^)n- 



(65) 
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and also at low temperatures: 

h.2\r.0 = ^,C- (66) 

We see that nonmagnetic disorder suppresses the orbital pair breaking and thus en- 
hances the upper critical field. 

In the general case, Hc2{T) can be calculated analytically only in the vicinity 
of the critical temperature using the Ginzburg-Landau free energy expansion. The 
results for the impurity response turn out to be nonuniversal, i.e. dependent on the 
pairing symmetry, the values of the intra- and interband coupling constants, and the 
densities of states in the helicity bands. li33J 



5 Spin susceptibility 

In this section we calculate the magnetic response of a noncentrosymmetric super- 
conductor, neglecting the orbital magnetic interaction and taking into account only 
the Zeeman coupling of the electron spins with a uniform external field H: 

where jis is the Bohr magneton. The components of the spin magnetic moment 
operator in the band representation have the following form: 

where ~ + 7^. 

The magnetization of the system is expressed in terms of the Green's functions 
as follows: J(i = tr'naA'(*)^A'V(*''^)' where ^(ifc,a)„) is the 

impurity-averaged 4x4 matrix Green's function in the presence of magnetic field 
(recall that the upper indices label the Nambu matrix components, see Sect. 12.1b . 
In a weak field, we have = Y^jXijHj, where is the spin susceptibility tensor. 
Treating the Zeeman interaction, Eq. (l67l l. as a small perturbation and expanding 
in powers of H, we obtain: 

y^^i = -r J:^; j:(TrM,(*:)^^(*:,*:';«„)M,(*:')^(*',*:;o^)),„,, (69) 
« it,*' 
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where Mj{k) = diag [mi{k),—mj{—k)]. The Green's functions here are unaveraged 
4x4 matrix Green's functions at zero field, and the trace is taken in both the 
electron-hole and helicity indices. 

In the clean case, one can evaluate expression (|69] l by summing over the Mat- 
subara frequencies first, followed by the momentum integration. The susceptibility 
tensor can be represented as Xij = xtj + Xij + Xij (Ref. EH), where 

are the intraband contributions, determined by the thermally-excited quasiparti- 
cles near the Fermi surfaces. Here Yx{k,T) = 2 Jq d^{—df /dEx) is the angle- 
resolved Yosida function, /(e) = (e^^^ + 1) ' is the Fermi-Dirac distribution func- 
tion. El {k) = ^J^^ + \Ax (^)p is the energy of quasiparticle excitations in the Ath 
band, and, as in the previous sections, the subscripts A in the Fermi-surface averages 
are omitted for brevity. The interband contribution is given by 

Since Xu is determined by all quasiparticles in the momentum-space shell "sand- 
wiched" between the Fermi surfaces, it is almost unchanged when the system un- 
dergoes the superconducting transition, in which only the electrons near the Fermi 
surface are affected. 

Collecting together the inter- and intraband contributions, we arrive at the fol- 
lowing expression for the spin susceptibility of a clean superconductor: 

Xij = Xij + IJ-INf Pa {rnVx ) ■ (72) 

At zero temperature there are no excitations (Yx = 0) and the intraband terms are 
absent, but the susceptibility still has a nonzero value given by Xij- The tempera- 
ture dependence of the susceptibility in the superconducting state at < T < is 
determined by the intraband terms, with the low-temperature behavior depending 
crucially on the magnitude of the SO coupling at the gap nodes. While in the fully 
gapped case the intraband susceptibility is exponentially small in all directions, in 
the presence of the lines of nodes it is proportional to either T or T^, depending 
on whether or not the zeros of Ax{k) coincide with those of y{k), see Ref. 1451 . 
In Fig.[3]the temperature dependence of Xij is plotted for a Rashba superconductor 
with y{k) = 7x(^ X z) and a cylindrical Fermi surface (referred to as the 2D model), 
and also for a cubic superconductor with y{k) = yok and a spherical Fermi surface 
(the 3D model). In both cases, the gaps in the two helicity bands are assumed to be 
isotropic and have the same magnitude. 
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Fig. 3 The clean-case temperature dependence of the transverse components of the susceptibihty 
for the 2D model (the solid line), and of all three components for the 3D model (the dashed line); 
Xp = IflgNp is the Pauli susceptibility. 

Now let us include the scalar disorder described by Eq. (|2]), or, equivalently, Eq. 
(|5]). After the impurity averaging, Eq. ( |69l ) is represented by a sum of the ladder di- 
agrams containing the average Green's functions ([TtI i. In contrast to the clean case, 
it is not possible to calculate the Matsubara sums before the momentum integrals. 
To make progress, one should add to and subtract from Eq. ( |69] l the normal-state 
susceptibility. 1251 It is easy to show that the latter is not affected by impurities and 
is therefore given by Xi^sj = Xij + I^b^f La Pi {fifj)^ see Eq. Then, 

Xij-XN.ij = -TY^' ^Y.{^TvM0Mj':^)i„,p- {TvM0NM0N),mp). (73) 
" kk! 

where ^^at is the unaveraged 4x4 matrix Green's function in the normal state with 
impurities, and the Matsubara summation is limited to the frequencies |a)„| < (Oc, at 
which the gap function is nonzero. 

Due to convergence of the expression on the right-hand side of Eq. ( TTTt , one can 
now do the momentum integrals first. The second term vanishes, while the calcu- 
lation of the ladder diagrams in the first term is facilitated by the observation that 
the integrals of the products of the Green's functions from different bands are small 
compared with their counterparts containing the Green's functions from the same 
band. The argument is similar to the one leading to Eq. dSTl i. For example, assuming 
that both the SO band splitting and the pairing are isotropic, i.e. ^±{k) = Eoik) ± 7 
and D± = D, we have 

jdSQGxG^X 1 rmax(a),.,r)l- 
max = max ^, — ^ ^ ^ <C 1 , 




where i2„ = ^ (b^ + In the same way, one can also obtain estimates for the 
momentum integrals containing anomalous Green's functions: 
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Eso 
max((B£.,F) 



max — . = max -;— — ^ , „^ ^ <c 1 , 

2 



~ — iiiaA. 



JdeoFxF_x 1 
max — - = max ttt^ 



Eso 



< 1. 



Thus we see that it is legitimate to keep only the same-band contributions to the 
impurity ladder on the right-hand side of Eq. ( |73] l. 

Following Ref. ||2T1 . we obtain the following expression for the spin susceptibil- 
ity: 

Z.=Z..,V + ^ri:^, (74) 
where X, are found from the equations 

- Y^i^uiXj +A2jjYj +Ai,jY*) = X(u, 



J 



i 

The notations here are as follows: 



(75) 



and Q„ ~ \J 6)}^ + \Dx P- Due to fast convergence, the Matsubara summation in Eq. 
(l74b can be extended to all frequencies. 



5.1 Residual susceptibility for isotropic pairing 

The general expression for the spin susceptibility, Eq. ( |74l l. is rather cumbersome. 
On the other hand, apphcation of our results to real noncentrosymmetric materials 
is complicated by the lack of information about the superconducting gap symmetry 
and the distribution of the pairing strength between the bands. As an illustration, we 
focus on the isotropic pairing model introduced in Sect. 13. II In this model, the order 
parameter magnitudes in the two bands are the same, but the relative phase can be 
either or n. While in the clean limit the spin susceptibility for both states is given 
by Eq. ( |72] |. the effects of impurities have to be analyzed separately. 

77+ = 77_ = T] . Solving equations (iTSl l we obtain, in the coordinate system in 
which (fifj) is diagonal, the following expression for the nonzero components of 
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the susceptibility tensor: 

MEl^l^l^f)nTy J'\ , ^ , (76) 

wherel]=(l-(f))r. 

We are particularly interested in the effect of disorder on the residual susceptibil- 
ity at zero temperature. In this limit, the Matsubara sum in Eq. ( l76T l can be replaced 
by a frequency integral, which gives 

Xp WoJ 
where tjq — {7z/e^)Tc is the gap magnitude at T = 0, and 



While the first two terms on the right-hand side of Eq. ( l77b represent the resid- 
ual susceptibility in the clean case, the last term describes the impurity effect. In 
a weakly-disordered superconductor, using the asymptotics <Pi{x) ~ nx/4, we find 
that the residual susceptibility increases linearly with disorder In the dirty limit, 
r ^ rjo, we have 'Pi{x) 1, therefore Xii{T = 0) approaches the normal-state 
value Xp- For the two simple band-structure models (2D and 3D) discussed earlier 
in this section, the Fermi-surface averages can be calculated analytically, and we 
obtain the results plotted in Fig.|4] 

Thus we see that, similarly to spin-orbit impurities in a usual centrosymmetric 
superconductor. 1461 scalar impurities in a noncentrosymmetric superconductor lead 
to an enhancement of the spin susceptibility at T = 0. Since the interband contri- 
bution is not sensitive to disorder, this effect can be attributed to an increase in the 
intraband susceptibilities. 

TJ+ = — ?]_ = 77. From Eqs. ( l74b and dTST l we obtain the nonzero components of 
the susceptibility tensor; 

^ = l-(f)7rrV— (78) 

Xp ^{co^ + ri2)y^-r iff 

We note that for a spherical 3D model with (ff) = 1/3 this expression has exactly 
the same form as the susceptibility of the superfluid ^He-B in aerogel, see Refs. Il40l 
andl|47l. 

At r = 0, the expression (iTST i takes the following form: 

^^^^^ = l-(1?) + (f)^2(0, (79) 



where 
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Fig. 4 The residual susceptibility at T = vs disorder strength for = r]_ . The solid line cor- 
responds to the transverse components in the 2D case {x~ = Xp and is disorder-independent), the 
dashed line - to the diagonal components in the 3D case. 




r/i7o 

Fig. 5 The residual susceptibility at T = vs disorder strength for ?]+ = — 7]_. The solid line 
corresponds to the transverse components in the 2D case (Xzz = Xp is disorder-independent), 
the dashed line - to the diagonal components in the 3D case. 



1- 



(/ + 1)3/2 



and y,nin = d{x— 1 ) \/ — 1 . The last term on the right-hand side of Eq. ( |79] l de- 
scribes the effect of impurities. According to Sect. 13.11 superconductivity is sup- 
pressed above the critical disorder strength = {'^ /'^e'^)TcO- For a given F, one 
should first obtain the gap magnitude from Eq. ( |37] ) and then calculate ^{x). In 
the weak disorder limit, we have 4>2(x) ~ (37r.x:/16)(l — {jj)), i.e. the residual sus- 
ceptibility increases linearly with disorder. At F i7, we have ^2{x) — > 1 and 
Xu{T = 0) Xp- The dependence of Xii{T = 0) on the disorder strength for the 2D 
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and 3D models is plotted in Fig.|5] As in the case t]+ — Tj-, the residual susceptibil- 
ity is enhanced by impurities. 



6 Conclusions 

Scalar disorder in noncentrosymmetric superconductors causes anisotropic mixing 
of the electron states in the bands split by the SO coupling. The critical temperature 
is generally suppressed by impurities, but this happens differently for conventional 
and unconventional pairing. For all types of unconventional pairing (which is de- 
fined as corresponding to a non-unity representation of the crystal point group, with 
vanishing Fermi-surface averages of the gap functions), the impurity effect on Tc is 
described by the universal Abrikosov-Gor'kov equation. The same is also true for 
certain types of conventional pairing, in particular the "protected" isotropic triplet 
state with 77+ = — 77_. Any deviation from the Abrikosov-Gor'kov curve, in particu- 
lar, an incomplete suppression of superconductivity by strong disorder, is a signature 
of conventional pairing symmetry. 

The impurity-induced mixing of singlet and triplet pairing channels makes the 
magnetic response of noncentrosymmetric superconductors with the SO coupling 
different from the centrosymmetric case. In an isotropic BCS-like model, the upper 
critical field Hc2 is enhanced by disorder at all temperatures, the magnitude of the 
effect depending on the SO coupling strength. In general, the effect of impurities on 
the slope of Hc2 is sensitive to the pairing symmetry and the band structure. 

Concerning the spin susceptibility, we found that scalar impurities in noncen- 
trosymmetric superconductors act similarly to spin-orbit impurities in centrosym- 
metric superconductors, in the sense that they enhance the residual susceptibility at 
T =0. The quantitative details again depend on the band structure, the anisotropy 
of the SO coupling, and the symmetry of the order parameter. 

Most of the experimental work on noncentrosymmetric superconductors has been 
done on CePt3Si. In this compound the Fermi surface is quite complicated and con- 
sists of multiple sheets. Il23l It is not known which of them are superconducting. The 
order parameter symmetry is not known either, although there is experimental evi- 
dence that there are lines of nodes in the gap. p8ll49ll50l The data on the impurity 
effects are controversial. The experimental samples seem to be rather clean, with 
the ratio of the elastic mean free path / to the coherence length ^0 ranging from 4 
(Ref. HI) to 10 — 27 (Ref. ISTl . There are indications that Tc is indeed suppressed 
by structural defects and/or impurities in some samples. ||50l On the other hand, the 
values of both the critical temperature and the upper critical field in polycrystalline 
samplesfl] are higher than in single crystals. [48 1 This is opposite to what has been 
observed in other unconventional superconductors and also disagrees with the theo- 
retical predictions, assuming that the polycrystals are intrinsically more disordered 
than the single crystals. In addition, the low-temperature behaviour of the penetra- 
tion depth in disordered samples is unusual ISOl and cannot be explained by existing 
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theoretical models. In order to resolve these issues, more systematic studies of the 
disorder effects in a wide range of impurity concentrations are needed. 
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